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ABSTRACT 

I examine the position of the ice line in circumbinary disks heated by steady mass accretion and 
stellar irradiation and compare with the critical semi-major axis, interior to which planetary orbits 
are unstable. There is a critical binary separation, dependent on the binary parameters and disk 
properties, for which the ice line lies within the critical semi-major axis for a given binary system. 
For an equal mass binary comprised of 1 M Q components, this critical separation is « 1.04 AU, and 
scales weakly with mass accretion rate and Rosscland mean opacity (oc [Mkr] 2 ' 9 ). Assuming a steady 
mass accretion rate of M ~ 10 -8 M Q yr _1 and a Rosseland mean opacity of kr, ~ 1 cm 2 g _1 , I show 
that > 80% of all binary systems with total masses M to t < 4.0 M Q have ice lines that lie interior 
to the critical semi-major axis. This suggests that rocky planets should not form in these systems, 
a prediction which can be tested by looking for planets around binaries with separations larger than 
the critical separation with Kepler (difficult) and with microlensing. 

Subject headings: planetary systems: protoplanctary disks - planetary systems: formation - stars: 
pre-main sequence - binaries: general 



1. INTRODUCTION 

The "ice line," or "snow line," is the point in proto- 
planctary disks at (and beyond) which the temperature 
is low enough that water-ice can condense from the neb- 
ular gas. Stellar irradiation, mass accretion, and general 
disk properties (e.g. opacity, surface mass density distri- 
bution) determine the position of the ice line. The posi- 
tion of the ice line is weakly dependent on the disk envi- 
ronment; stellar irradiation, mass accretion, and general 
disk properties determine the position. The temperature 
threshold below which ice can form is ~ 145 — 170 K, 
depen ding on the partial pressure of the nebular water 
vapor (|Podolak k Zuckerll200l iLecar et al J 120061) . The 
ice line is important in the context of planetary forma- 
tion because it separates the inner region of rocky planet 
formation from the outer region of gaseous and icy plan et 
formation (jlda k Linll20ol [Kennedy k Kenvonl 12001 ) . 

Gas giant formation by core accretion occurs when pro- 
toplanet co re masses are suffic ient to attract gas from 
the nebula. Ilkoma et al.l (|2000l) show that the timescale 
for gas giant format ion is set by the core mass (see also 
IHubickvi et al.|[2005h . Thus, if the core mass is too small, 
it will not be able to accrete enough gas to form a gas gi- 
ant before the gas is removed. Past the ice line, condensa- 
tion of water ices increase the surface density of the disk 
by a factor of ~ 3, increasing the isolation (core) mass by 
a factor of ~ 5 (Kennedy & Kcnyon 2008). This increase 
in surface density beyond the ice line allows for gas giant 
formation by core accretion within the short lifetime of 
disk gas , believed to be somewhere in the range of ~1- 
10 Myr ijZuckerman et alJll995HPascucci et al.)l2006D . It 
is therefore important to know where the snow line lies 
in protodisks, as it holds much information about where 
gas giants are formed and how their orbits evolve in time. 

Only recently have circumbinary planets been discov- 
ered. There are currently six known circumbinary plan- 
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ets in five different systems found by Kepler (jDovle et al.l 
l20TllWelsh et^l20rllOrosz et al.N2012allbh . Addition- 
ally, there are 12 eclipsing, close-compact binary systems 
with white dwarf or hot sdB primaries showing appar- 
ent period variations that could perhaps be explained by 
circumbinary giant planets. These suspected giant plan- 
ets were probably not formed primordially, but rather 
it is more likely the y formed later, after commo n en- 
velope evolution (see lZorotovic k Schreiberl (|2013j ) for a 
detailed discussion of these systems). Nevertheless, the 
unambiguous discovery of circumbinary planets by Ke- 
pler has motivated a growing body of literature regarding 
planetary formation in circumbinary disks, including this 
study. 

In this paper, I study ice line positions in circumbinary 
disks and investigate the parameters that determine their 
locations. In § [2] I describe a simple circumbinary disk 
model, and in § [3] I calculate ice line positions and de- 
scribe results. Conclusions are presented in § 0J 



2. MODEL 



2.1. 



Disk Model and Stellar Parameters 

Similar to the treatment of ILecar et al.l (|2006f ) in the 
case of a single star, I consider a circumbinary disk 
heated by steady mass accretion as well as stellar irradi- 
ation from the binary. I adopt a cir cumbinary disk with 
a flare d geometry, as described in iChiang k Goldreichl 
(|1997f) . The surface mass density distribution is S = 

So (a /AUr 3/2 , with S = 10 3 g cm" 2 (|Weidenschillind 
[197% . 

Several simulations have shown that the inner edge 
of circumbinary disks should be truncated between 
1.8 o,bin and 2.6 Obin, where dbin is the binary semi- 
major axis, due to tidal torques exerted by the 
binar y (lArtvmowicz et all 119911: [Artvmowicz k Lubowl 



1991 iG^mther fc Kiev! 120021: IPierens k Nelson! I2007t 



MacFadven k Milosavlievid I2008TT There has been 
general agreement with these predictions from ob- 
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servations of directly imaged circumbinary disks, 
most notabl y tha t found a round GG Tauri by 
Dutrev et all (119941) (see also iBeust fe Dutrevl 120051: 



Duvert et al.lll998tllreland fc Kra us 2008). Furthermore, 
Holman fc Wiegertl (|1999T ) have shown that there is a 
critical semi-major axis for planets orbiting a binary. 
Within this critical semi-major axis, the planets are too 
strongl y affected by the binary a nd their orbits are un- 
stable. iHolman fe Wiegertl (|1999f ) provide the relation 

a c = [(1.60 ± 0.04) + (5.10 ± 0.05) e 

+ (-2.22 ± 0.11) e 2 + (4.12 ± 0.09) fj, 

+ (-4.27±0.17) e/j + (-5.09 ±0.11) fx 2 

+ (4.61 ±0.36) eV] abi„ , (1) 

where a c is the critical semi-major axis, (i = 
va-ij (mi + 77i2 ) is the mass ratio of the binary, and e is 
the eccentricity of the binary. For an equal mass binary 
on a circular orbit (/i = 1/2, e = 0), this relation yields 
a c ~ 2.4 abin, consistent with the truncation radius ex- 
pected for circumbinary disks. 

I construct pre-main sequence binary systems using 
stella r masses, t e mper atures, and radii from the mod- 
els of Sies s et al.1 ()2000l ) at an age of one-tenth the stars' 
zero-age main sequence (ZAMS), tzAMs/10- I examine 
both equal and unequal mass binaries, with component 
masses ranging from 0.5 — 5.0 Mq. 

2.2. Midplane Temperature Profiles 

Accretion Heating: As in lLecar et al.l (|2006l ). I adopt 
the lLvnden-Bell fc Pringlel (|1974h prescription for steady 
accretion 
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To generalize this to the case of a binary, I replace M* 
with the total mass of the binary, M tot . Note that I 
ignore the \/ R*/r term usually associated with equation 
[2 as it is a very small correction for this study due to 
the truncated inner edges of circumbinary disks. I adopt 
a constant value of Af = 10 -8 M Q yr -1 for the mass 
accretion ratcQ. 

The temperature due to accretional flux diffusing in an 
optically thick medium from the midplane to the surface 
of the disk is 
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where tr is the Rosseland optical depth 
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which I approximate as tr ~ kr£ (r) /2, where kr is the 
Rosseland mean opa city. I adopt the opa city value of 
kr ps 1 cm 2 g _1 from lD'Alessio etahl (|2001l ). using their 
dust model with a maximum grain size of a max = 1 mm 
and grain size distribution power law slope of p = 3.5. 



1 IHartmann et all 11199 9) found a median accretion rate for T 
Tauri stars of ~ 10 — 8 Mq yr -1 , with a spread as large as one 
order of magnitude, in the Taurus and Chamaelcon I molecular 
cloud complexes. I discuss the sensitivity of my results to this 
choice in § [3] 



Uang-Condell fc Sasselovl (|2004l ) showed that the depen- 
dence of this opacity on the temperature structure of a 
disk is small between 100 — 300 K, so I assume that kr 
is constant throughout the regions of the disk I examine. 

Irradiation: Around a binary system the ice line 
will oscillate due to the orbital motion of the irradi- 
ating sources. The extrema bounding the ice line de- 
pends on the stellar temperatures (T^^T*^), stellar 
radii (R*,i, -R*,2)i stellar masses (M^i, M* j2 ), and the ra- 
tio of the orbital period of the binary to the cooling time 
for the disk material (r = -PbmAcooi)- Figure [TJ shows 
instantaneous orbital configurations of the binary that 
produce the two extrema between which the ice line will 
lie. The semi-major axes of the primary and secondary 
I define to be ab,i and <2b,2, respectively. 

For systems where the cooling time of the disk material 
is much longer than the binary orbital period (r <C 1), 
the ice line will be static, with a position correspond- 
ing to the maximum ice line position of the r > 1 sce- 
nario. Since I am interested in the ice line as it relates 
to planet formation, I only consider the maximum po- 
sition of the ice line in this study becase in the r > 1 
case, planetesimals that form from condensed water-ice 
within this maximum ice line will be reheated as the bi- 
nary orbits. Since binary orbital timescales are much 
shorter than timescales for planet formation, water-ice- 
based planetesimals within the maximum ice line will be 
too short lived to form planets. 

For configuration A shown in Figurc[U the flux incident 
at a point P in the disk is 
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the irradiation strikes the disk, given by 
0.4i?* 8 ( T, 
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For the parameters adopted, tr > 1 for both the in- 
coming radiation and the radiation from the midplane 
for a < 63 AU. In this region, the inner disk tem- 
perature of a passive disk is Tj rri A = J 1 e /2 1 ^ 4 , where 

T e = ( Fa/ct) 1 ^ 4 , which is indepe ndent of disk surface 
density (jChiang fc Goldreicb][r997l ). 

For configuration B, pictured in Figure [TJ I assume 
that the primary fully eclipses the secondary, such that 
the only flux the disk receives at the point P is from the 
primary. The flux incident on the disk is 



Fn 



(t)U 



- a b ,i (1 ± e) 



(7) 



and the temperature of the interior of the disk is Ti ri i B = 

(Fb/2<7) 1 ^ 4 . Configuration C, similar in geometry to B, 
shown in Figure [TJ I assume the secondary fully eclipses 
the primary, yielding an internal disk temperature of 
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Config A Configs B/C 




Inner Disk Edge Inner Disk Edge 



Figure 1. Cartoon showing the various geometric configurations of the binary, each leading to a different position of the ice line. The sign 
of the eccentricity terms depends on whether the stars are at apastron (+) or periastron ( — ). The star labeled as "1" is the primary and 
the star labeled as "2" is the secondary. Configurations B and C are geometrically similar, just with the primary and secondary switched. 
I assume that in configurations B and C, the star nearest to the point P in the circumbinary disk fully eclipses its companion. 

a/4 



7irr,c = (-Fc/2er) 1/ ' 1 , where Fc is given by equation (0 
but with i?*^ — > i?*.2, T*,i — > r* 2, a\ — > ait and 
Ob.i ~> ib, 2- Note that in the case of equal mass binaries, 
configurations B and C are degenerate. 

Combining heating due to accretion and irradiation, 
the midplane temperature of this circumbinary disk 
model is given by 



mid 
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3. RESULTS 



I find the position of the ice line by determining the 
point in the circumbinary disk where eq uation ([8| is 
equal to 160 K. Originally, [HayashJ ( |198lD took the ice 
line condensation tempe rature to be 170 K, however, 
iPodolak fc Zuckeil (|2004[ ) showed that the ice line tem- 
perature can be as low as ~ 145 K. Rather than bal- 
ancing the evaporative cooling and condensation heating 
for the grains to find the position of the ice line, I adopt 
an ice line temperature of 160 K, since the position of 
the ice line will not change significantly for condensation 
temperatures in this range. 

Figure [2] shows the maximum ice line position relative 
to the critical semi-major axis, interior to which orbits 
are unstable (see equation (fl) ) , as a function of the bi- 
nary separation for both equal mass binaries (top panel) 
and unequal mass binaries (bottom panel) on circular or- 
bits. The unequal mass systems are constructed from a 
0.5 M star as a secondary. For the /i ~ 0.2, \i ~ 0.25, 
\i ~ 0.33, and (i ~ 0.42 systems, the primary has a mass 
of 2 M Q , 1.5 M Q , 1 M Q , and 0.7 M Q , respectively. 

In all systems I consider that have component masses 
< 2.0 Mq, heating by accretion dominates the deter- 
mination of the ice line position. Assuming that for 
separations from the ice line inwards there is only ac- 
cretional heating, equation ([5]) can be approximated as 
Tiaid = T acc . Combining this with equations ([5]) and J3]), 
along with the fact that tr 3> 2/3 in this regime, 
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tal binary mass Af to t is independent of the binary sep- 
aration. Normalizing this ice line position to the criti- 
cal semi-major axis, which is proportional to the binary 
separation for a given mass ratio and eccentricity, it is 
clear that a- lcc /a c oc a^ in - This proportionality holds for 
any binary eccentricity and mass ratio, provided that no 
component of the binary has a mass > 2.0 Mq, so that 
accretion is the dominant heating mechanism out to the 
position of the ice line. 



2.0 



o 



1.5- 



o 




0.5 M s 
0.7 M 
1.0 M s 
1.5 M 
2.0 M s 
2.5 M 
3.0 M s 
4.0 M s 
5.0 M s 



0.20 (M,=2.0 M G 
0.25 (M,=1.5 M e 
0.33 M,=1.0 M G 
■ 0.42 M,=0.7 WU 



=0.5 M £ 
=0.5 M a 
=0.5 M £ 
=0.5 M ff 




which shows the position of the ice line for a given to- 
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Figure 2. Position of the ice line relative to the critical semi- 
major axis as a function of binary separation for binaries on cir- 
cular orbits. The top panel shows curves for equal mass bina- 
ries ()i = 1/2), with the mass of each component labeled in the 
legend, while the bottom panel allows for various mass ratios, 
li = Mil (Mi + M 2 ). 
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The most important information to be gleaned from 
Figure [2] is that in each case, there is a critical binary 
separation for which a; cc /a c = 1, which I define as a"^ = 
Cfbin (flicc/cJc = 1)- For binary separations equal to and 
larger than the ice line is always interior to the 

inner disk edge. Analytically, employing equation ([9]) 
and the expression for the critical semi-major axis as a 
function of binary separation for equal mass binaries on 
circular orbits (a c ~ 2.4at,i n ), 
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(10) 



which is good to within a few percent for all the equal 
mass systems shown in Figure [5] with component masses 
< 2.0 M Q . Setting af c c e c /a c = 1 yields an expression for 

,.crit 
"bin- 

The equal mass binaries with component masses > 
2.5 Mq heat their disks primarily by stellar irradiation. 
Assuming T m id = T m = Xi CC for equal mass binaries on 
circular orbits in configuration A shown in Figure Q] yields 
no closed- form solution for a"* c /a c . However, an analyt- 
ical expression for ajjJJJ resulting from heating by stellar 
irradiation can be obtained: 
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If the binaries arc allowed to have eccentric orbits, the 
critical semi-major axis, a c , grows via equation ([1]). All 
unequal mass binaries considered and the equal mass bi- 
naries with component masses < 1.0 Mq remain in the 
accretion-dominated regime for all eccentricities. How- 
ever, for equal mass binaries with component masses 
> 1.0 M Q , irradiation dominates at distances near the 
ice line for binary separations equal to the critical value, 
a bin- l n both the accretion- and irradiation-dominated 
regimes, since a c always increases, a£"* decreases with 
eccentricity. Figure [3] shows the change in a"^ as a func- 
tion of binary eccentricity for systems of both equal and 
unequal mass binaries. 

For the systems in the accretion-dominated regime, the 
change in a£"* as a function of binary eccentricity can be 
expressed analytically by combining equations ([9]) and 
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where 4> (e, £t) is the coefficient for the critical semi-major 
axis, given by the expression in the brackets of equation 
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Figure 3. Change in the critical binary separation, aKjJ, defined 
to be the binary separation where the ice line is equal to or interior 
to the critical semi-major axis, a lcc < a c , relative to circular orbits 
as a function of binary eccentricity. The top panel shows that 
for equal mass binaries and the bottom panel shows the same for 
unequal mass binaries. The diamonds in each plot are calculated 
points at steps of Ae = 0.1, color coded by component mass or 
mass ratio. The systems in the accretion-do min ated regime have 
analytical lines drawn according to equation 1121 . 

|T]) . For the accretion-dominated systems in Figure El I 
plot these analytic curves. I find no simple, analytic form 
for Aa£j^ in the irradiation-dominated regime. The main 
point from Figure [3] is that binary eccentricity serves to 
decrease a™£ m au cases. 

I also examine how the maximum ice line changes in 
time, as the binary components contract to the main 
sequence, for equal mass binaries o n circular orbits. U s- 
ing the stellar evolution models of iSiess et al.l (|2000f ) . I 
calculated ag^ for each system as a function of time, as 
the binary components approach the main sequence. For 
equal mass binaries consisting of stars less massive than 
~ 2 M Q , the temperature and luminosity do not increase 
enough to dominate disk heating over steady mass accre- 
tion until they reach ZAMS. However, even at ZAMS, the 
ice line is not increased much from the initial position at 
formation time. Furthermore, by the time stars less mas- 
sive than ~ 2 M Q reach the ZAMS, their disks will have 
dissipated and any planet formation should have already 
ocurred. Conversely, equal mass binaries consisting of 
stars more massive than ~ 2 M Q will indeed have their 
ice line locations significantly increased as they contract 
towards the main sequence. The effective temperature of 
such massive stars increases rapidly as they approach the 
main sequence, and stellar irradiation quickly becomes 
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the dominant heating mechanism for their disks. Since 
izAMS for massive stars is much less than disk lifetimes, 
any planet formation around massive star binaries will 
necessarily continue after they reach the main sequence. 
For the purposes of this study, I thus restrict my conclu- 
sions to binaries consisting of low to intermediate mass 
components, M* < 2 M Q . 

It is also worth noting that the mass accretion 
rate among T Tauri stars declines as the stars age 
(|Hartmann et al.l [ l9 98). However, since af c c e c oc M 2 / 9 , 
a declining mass accretion rate will serve to move the 
ice line inwards. iLecar et al.l (|2006t ) pointed out that 
the disk lifetime scales as 1/M and that the lifetime 
of the minimum-mass solar nebula disk for a constant 
M = 1CT 8 Mq yr _1 is 2 x 10 6 yr, or roughly the life- 
time of the nebular gas. Accretion rates higher than 
M Q yr -1 would further constrain the (already dif- 
ficult to explain) timescales for giant planet formation. 
Nonetheless, the scaling of the ice line in the accretion- 
dominated regime with the mass accretion rate is weak 
and thus variations of an order of magnitude would not 
significantly change the results of this study. 

It is interesting that there is a maximum binary separa- 
tion (i.e. (Zbin) f° r w hi cn systems with separations larger 
than this maximum value have ice lines that lie interior to 
the critical semi-major axis. Circumbinary disks around 
systems with a D m > a£in w ^ therefore have condensed 
water-ice throughout the entire disk, meaning that rocky 
planets should not form in such systems. Using the prop- 
erties of the observed binary population, I calculate the 
fr action of these b i nary s ystems with a D m > a bin - 

IRaghavan et al.l (|2010D determine the orbital period 
distribution of solar-type stars with stellar companions 
to be consistent with a Gaussian, with a mean of log P = 
5.03 and eri og p = 2.28, where P i s in days. I assume that 
this period distribution found by IRaghavan et al.l (|2010l ) 
is the same for all equal mass binaries with component 
masses in the range I consider, 0.5 M Q < M* < 2.0 M Q . 
I then compute the the periods corresponding to a b "' for 
equal mass binaries with component masses in this range 
and integrate over a Gaussian distribution with the above 
parameters to estimate the fraction of binaries with ice 
lines interior to their inner disk edge. 

For equal mass binaries with total masses M to t < 
4.0 M , the percentage of binary systems that have ice 
lines interior to the critical semi-major axis, at izAMs/10 
and later, ranges from 81% (M tot = 4 M Q ) to 84% 
(M to t = 1 Mq). These fractions represent lower limits 
because I conservatively choose the "maximum" ice line 
in each system. Furthermore, a fraction of 80% is a con- 
servative lower limit for unequal mass binaries since ajj-jj 
for the lowest mass ratio system I consider (/x ~ 0.2) is 
approximately the same as that of the equal mass binary 
with component masses of 2.0 M Q . 

4. CONCLUSION 

I have shown that ice lines lie interior to the critical 
semi-major axis for > 80% of all equal and unequal mass 
binaries with components less massive than M* < 2 M Q . 
This suggests that rocky planets should not form in these 



systems; only gas giant and icy planets. 

These results can be tested by finding more circumbi- 
nary planets, specifically planets around binaries with 
separations larger than the critical separation. This is 
possible, albeit difficult, with Kepler, because there are 
fewer eclipsing binaries at larger separations and plan- 
etary transit probabilities decrease with orbital period. 
Microlensing can also detect circumbinary planets at a 
range of binary separations and seems to be the most 
promising method for testing the results of this study 
(|Hai][2T)08h . 

I gratefully acknowledge helpful conversations with 
Todd Thompson and Scott Gaudi about this project, as 
well as critical readings of the text. 
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